Abstract. Using the standard infall model of Galactic 
chemical evolution, we explore the origin of carbon and 
calculate the abundance evolution of CNO elements for 
8 different models of stellar nucleosynthesis yields. The 
results show that, in the early stage of the Galaxy, mas- 
sive stars are the main producer of carbon, and that as 
our Galaxy evolves to the late stage, the longer lived 
intermediate- and low-mass stars play an increasingly im- 
portant role, while at the same time, metal-rich Wolf- 
Rayet stars eject a significant amount of carbon into the 
ISM by radiative-driven stellar winds. However, from the 
< T^) • present published nucleosynthesis yields we cannot distin- 
guish whether the main source of carbon in the late Galac- 
tic stage is just the massive stars (M > 8M ) alone, or 
just the intermediate-, low-mass stars and M < 40-M© 
massive stars that do not go through the Wolf-Rayet 
stage. The 12 C(a,7) 16 reaction rate is very important 
in the stellar nucleosynthesis calculations: a lower rate 
will give a higher yield of carbon. The contribution to ni- 
trogen is dominated by intermediate- and low-mass stars, 
and the secondary source of massive stars cannot explain 
the observed [N/Fe] in metal-poor stars. Most of oxygen 
is produced by massive stars. The fact that a higher O 
abundance in metal-poor stars is derived from the O I 
7771-7775 A triplet than from the forbidden [Oi] line 
at 6300 A poses a problem. 
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1. Introduction 

The abundance trends of the CNO elements are crucial to 
study the Galactic chemical evolution (hereafter GCE). 
These elements play important roles in the stellar interior 
as sources of opacity and energy production through the 
CNO cycle, and thus affect the star's lifetime, its position 
in the HR diagram, and its heavy-element yields. Carbon 
is especially important because it is the first element that 
is synthesized out of the primordial elements (H, He, Li). 
The main site of carbon synthesis has been a subject of 



the Galactic disk (—1.6 < [Fe/H] < +0.2) and some ir- 
regular galaxies. They suggested that carbon enrichment 
is caused mainly by superwinds from metal-rich, massive 



argument for many years. Burbidge et al. (1957) suggested 
that carbon is provided by m ass lo ss from red giants and 
supergiants. Iben & Truran ( 1978 ) concluded from ther- 
mally pulsing models that intermediate-mass and high- 
mass stars contribute carbon in roughly equal amounts. 
More recent theoretical and observational studies include 
the following. 



Timmes et al. ( 1995 ) used the carbon yields of Type 
II supernov a (her eafter SNII) explosion from Woosley 
& Weaver ( 1995 ) (hereafter WW) and the yields of 
intermediate -, low -mass stars (hereafter ILMS) from Ren- 
zini & Voli ( 1981 ) (hereafter RV) to calculate the abun- 
dance evolution of carbon. They concluded that the ILMS 
contribute a significant amount of carbon to the interstel- 
lar medium (hereafter ISM). Chiappini et al. ( |1997[ ) ob- 
tained similar results with their two-infall model. Ober- 
hummer et al. ( |2000| ) reviewed the ILMS as the major 
producer of carbon. Kobulnicky & Skillman (1998) mea- 
sured the C/O and N/O ratios of three metal-poor galax- 
ies with similar metallicities using HST ultraviolet and 
ground-based optical spectroscopy. Their results implied 
that the mechanisms of C and N production are coupled, 
and most of both C and N in galaxies originates in low- 
to intermediate-mass stars. 



Prantzos et al. (1994) suggested that if the duration 
of the halo phase is as long as ~ 1 — 2 Gyr, then massive 
stars with wind-driven mass loss could be the main car- 
bon source during the whole galactic evolution. Gustafs- 
son et al. (1999) adopted an analytical GCE model to 
calculate and compare the carbon abundance evolution in 
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stars. Karlsson et al. (1999) reported similar results. How- 
ever, Gustafsson & Ryde ( 2000) ) wrote that the source of 
carbon is still not clear. Henry et al. (2000) analyzed the 
behavior of the N/O and C/O abundance ratios as func- 
tions of mctallicity in the Galactic and extragalactic Hn 
regions. Their results showed that C and N originate from 
separate production sites and are decoupled from one an- 
other; and they suggested that massive stars (M > 8Mq) 
dominate the production of carbon, while intermediate- 
mass stars between 4 and 8M Q dominate the nitrogen pro- 
duction. Carigi ( [2000 ) adopted different sets of yields to 
calculate the C/O vs. Galactic age relation, and suggested 
massive stars to be the main contributor of carbon. 

From calculated abundance gradients, Hou et al. 
( 2000D argued that carbon comes from stellar winds of 
massive stars and ILMS. Goswami & Prantzos (2000) us- 
ing their new GCE model, concluded that the yields in 
WW are not enough to explain the observed carbon abun- 
dance, and that some other sources, Wolf-Rayet (hereafter 
W-R) stars or ILMS, are needed to explain the observed 
[C/Fe] vs. [Fe/H] curve. Garnett et al. ( 1999| ) measured 
the gas-phase abundance ratio C/O in six Hn regions in 
the spiral galaxies M 101 and NGC 2403, based on ultravi- 
olet spectroscopy from the HST. These authors suggested 
that stellar winds in massive stars have an important ef- 
fect on carbon and oxygen, and their analysis showed that 
the source of carbon depended on the assumed UV redden- 
ing function. When a flatter UV reddening (Ry — 5) was 
used, the data were consistent with the model that used 
the yields of massive stars from WW, while with a steeper 
value (i?y = 3) the data became more co nsiste nt with the 
model that used the yields of Maeder (1992) (hereafter 
M92). 

In summary, it is still not clear whether carbon comes 
from the SNII explosions of massive stars, from stellar 
winds of W-R stars or from ILMS. In this paper we try to 
address this question from two aspects: (1) by using GCE 
model calculations with different sets of nucleosynthesis 
yields; (2) by using the observations of the large sample 
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of dwarfs in our Galaxy as one of the main bservational 
constraints. 

For the first aspect, GCE modelling is an excellent tool 
to study the evolution of elements of the whole Galaxy 
and other galaxies. The main parameters for setting up 
such a model include the form of the infall, the star for- 
mation rate (SFR), the initial mass function (IMF) and 
the nucleosynthesis yields. The exponential form of infall 
(or Gaussian infall form: Chang et al. 1999| and Prant- 
zos & Silk 1998) has been found to be able to solve the 
G-dwarfs problem very well. The SFR and IMF can be 
determined fro m the observed present-day mass function 
(PDMF) (Pagel |l997| ). The most uncertain factor of GCE 
model is the nucleosynthesis yields, which largely decide 
the evolutionary trends of the elements. Different research 
groups have calculated the yields of stellar nucleosynthe- 
sis, but the results are rather different due to the use of 
different input parameters. And different yields will result 
in different evolution trends. 

We set up a simple GCE model, like the one of Timmes 
et al. (1995), which fits some important observations. We 



choose the following different nucleosynthesis yields pub- 
lished recently in our calculations: for the RV yields of 
intermediate- and low-mass single stars, van den Hoek 
& Groenewegen (1997) (hereafter VG), Marigo et al. 



fll996Q (hereafter MBC1), Marigo et al. ( |1998| ) (hereafter 
MBC2) (hereafter, MBC repre sents for combined MBC1 
and MBC2) and Marigo ( p00l| ) (hereafter M2K); for the 
WW yields of massive single stars, Nomoto et al. (1997a) 
(hereafter N97), M92 and Portinari et al. ( |1998| ) (hereafter 
PCB); for the yields of binaries through SNIa explosion, 
Nomoto et al. ( |1997b ). Then we compared and analyzed 
the different evolutionary behavior of the CNO elements, 
particularly C, calculated on the basis of these different 
sets of yields. 

In the early stage of our Galaxy, massive stars are the 
main producers of carbon. As our Galaxy evolves to a late 
stage, ILMS play an increasingly important role in the en- 
richment of the ISM due to their longer lifetimes. At the 
same time, the high metallicity W-R stars eject significant 
amounts of carbon into the ISM by radiative-driven stellar 
wind (the metallicity influences the opacities in the outer 
layers of massive stars, hence their mass loss rate through 
winds). Which is the main producer of carbon in the late 
stage of our Galaxy? We attempt to answer in the follow- 
ing way. We divide the massive stars into two mass ranges, 
M > 4OM and M < 40M©. And assume that the former 
does, and the latter does not, go through the W-R stage. 
We then compare the contributions by the W-R stars and 
ILMS to carbon in the late Galactic stage. 

In the second part of our work, we pay particu- 
lar attention to the observations of field dwarfs in our 
Galaxy. Their atmospheres are assumed to represent per- 
fect records of the chemical compositions of the ISM at the 
time of their formation, and we have available the abun- 
dances for a large sample of field dwarfs in the solar neigh- 



borhood. In addition, the metallicity distribution of G 
dwarfs is very important, because the G dwarfs have life- 
times as long as the Galaxy, their metallicity distribution 
reflects the evolution of our Galaxy. The age-metallicity 
relation is significant, too, for reflecting the Galactic evo- 
lution. 

This paper is organized as follows. In Sect. 2, we an- 
alyze the observational data. The GCE model and the 
relevant parameters are presented in Sect. 3. Different nu- 
cleosynthesis yields are analyzed in Sect. 4. Results of the 
analysis are given in Sect. 5—7. The main conclusions are 
summarized in Sect. 8. 

2. Observational data 

Observations of dwarf stars in the general field are very 
important to the evolution of chemical abundance as their 
atmospheres are thought to represent perfect records of 
the chemical compositions of the ISM at the time they 
formed. 

2.1. [C/Fe] vs. [Fe/Hj 

Carbon abundance can be calculated from the CH lines 
in the blue and ultraviolet (UV), the infrared CO emis- 
sion, and a few red and near-infrared C I lines and [C i] 
forbidden lines. The [Cl] line of low excitation poten- 
tial being not so tightly dependent on NLTE considera- 
tions, the values derived from it are less dependent on the 
adopted stellar effective temperature than those from the 
high-excitation C I lines and CH lines. However, the small 
equivalent width (EW) of the [C i] line makes the measure- 
ment difficult. And the [Cl] line at A 8727 A is supposed 
to be blended with a weak Fel line. Both the [Cl] and 
C I lines become unmeasurably weak in metal-poor stars 
([Fe/H] < -1). 

The CH line in the UV can be used to obtain the 



C abundance of halo and disk stars. Laird ( |1985| ) pub- 
lished an extensive study of the C abundances of 116 field 
dwarfs, 10 faint field giants and 3 Hyades dwarfs in the 
range -2.45 < [Fe/H] < +0.5. They suggested that the 
[C/Fe] were solar in the whole [Fe/H] range of the sam- 
ples. Tomkin et al. ( 1986| ) analyzed the spectra of 32 halo 
dwarfs, and derived [C/Fe] rj —0.2, independent of metal- 
licity for stars with [Fe/H] > —2. And their data showed 
a rise of C abundance in more metal-poor stars: [C/Fe] 
w +0.2 at [Fe/H] w -2.5. Carbon et al. ( |1987| ) derived the 
C abundances of 83 halo dwarfs, which showed that [C/Fe] 
is constant over the range of —2.5 < [Fe/H] < —0.75, and 
they found an uptur n at the very low metallicity value. 
Wheeler et al. ( 1989 ) reviewed these observed trends of 
[C/Fe] vs. [Fe/H]. 

Tomkin et al. (1992) derived the C abundances of 



halo and disk stars using the Cl A 9100 A line. Their re- 
sults showed that the C abundances are +0.2 higher than 
those derived from the CH lines. Andersson & Edvards- 
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son Q1994D (using [Ci] 8727 A), Tomkin et al . ( |1995[ ) (us- 
ing C I 7100 A) and Gustafsson et al. ( |1999| ) (using [C i] 
8727 A) obtained the C abundances for large samples of 
disk dwarfs (comprising 85, 105 and 80 stars respectively, 
with some in common). Their results showed that [C/Fe] 
decreased with increasing [Fc/H] in the Galactic disk re- 
gion. Using CH and C2 lines, Carretta et al. (2000) an- 
alyzed a sample of 19 stars, and reanalyzed the stars of 
Tomkin et al. Ql992j ) and Edvardsson et al. ( |1993| ). Their 
results showed that [C/Fe] is roughly solar in the whole 
metallicity range of [Fe/H] > —2.5. 

Considering the large scatter, we choose to take all the 
above-mentioned observations of C abundance as observa- 
tional constraints on our GCE calculations. 



2.2. [O/Fe] vs. [Fe/H] 

Oxygen is the third most abundant element in the uni- 
verse, and it is very important for the study of GCE. 
In general, oxygen is produced by massive stars through 
SNII explosions, like other a elements. Thus, the obser- 
vations of oxygen, especially in metal-poor region, give a 
measure of the contribution by massive stars. Some pa- 
rameters required in GCE calculations can be determined 



by the observed [O/Fe] in the stars (Prantzos et al. 1994). 

The oxygen abundance of stars can be derived from the 
Oi AA 7771-7775 A triplet, the [Oi] A 6300 A forbidden 
line and OH lines in the UV region. In main-sequence, 
metal-poor stars, the [0 1] lines are extremely weak and 
are difficult to use for this purpose. Nonetheless, sensi- 
tive detectors, high-resolution, and long exposures have 
enabled several teams (Spite & Spite 1991; Nissen & Ed- 



vardsson 1992) to measure the EWs of some dwarfs with 
[Fe/H] < —0.1. However, the [Oi] line becomes too weak 
to measure at [Fe/H] < —1.3, and to provide useful in- 
formation on the trend of [O/Fe] vs. [Fe/H]. Thus, for 
the metal-poor stars, the 1 triplet lines are the preferred 
choice in dwarfs because they are of measurable strength 
for a wide range of stellar temperatures and they lie in 
a fairly clear portion of the stellar spectrum. Molecular 
OH lines in the near-ultraviolet at AA 3060-3200 A, are 
observed in both giants and dwarfs with T c g < 6500 K 
(OH lines are destroyed at higher temperatures), and can 
be used to derive the O abundance. 

An interesting result is that there is a systematic dif- 
ference of about 0.5 dex between the O abundances de- 
rived from the permitted and forbidden lines. The [O 1] 
6300 line provides about 0.4—0.5 dex of the lower [O/Fc] 
values in metal-poor stars (e.g. Spite & Spite 1991, Spies- 
man & Wallerstein |1991| , Fulbright & Kraf t |1999| and 



Nissen et al. 2000 for dwarfs; Barbuy 



Erdelyi-Mendes 1989 and Sneden et al 



1988 



1991 



Barbuy & 
for giants). 

The 1 triplet and OH lines, on the other hand, produce 
the higher [O/F e] valu es, re aching +1.0 at [Fe/ H] ss —3 
(Abia & Reb olo |l989| ; King |l994j Ni ssen et al. |l994| ; Is- 
raelian et al. 1998] ; Boesgaard et al. 1999] ; Mishenina et 



al. 2000). Perhaps the higher O abundance obtained from 
the 1 triplet was caused by some NLTE effect (Mishenina 
et al. 200C and references therein). Kiselman ( |1991| ) has 
performed a NLTE analysis of Oi A 7774 A and found 
significant NLTE corrections, up to 0.4 dex. However, 



Takeda (|1994|) found the NLTE effects to be small in 
metal-poor stars. Another source of discrepancy may be 
the fact that the O/Fe ratio is very sensitive to T e g (Nis- 
sen & Edvardsson 1992j ). King (1992) determined a re- 
duction in the temperature of subgiants of 150-200 K, and 
obtained similar O abundances from the A 7774 A and 



A 6300 A lines. King & Boesgaard ( 1995 ) obtained the O 
abundances from both the [0 1] A 6300 A and 1 A 7774 
A lines for a sample of metal-rich F and G dwarfs. They 
found that for T c g < 6200 — 6300 K, no systematic differ- 
ence exists in the abundances determined from the AA 6300 
and 7774 A lines. For T cB > 6200 - 6300 K, however, 
the A 7774 A abundance is substantially higher than the 
A 6300 A abundance. 

Thus, many papers have demonstrated that there is a 
systematic difference between the values of O abundance 
derived from the forbidden and permitted lines. Possible 
reasons are as follows. The excitation potential of the 1 
triplet is very high (9.15 eV), the line is formed deep down 
in the atmosphere, and the abundance determination has 
a large error; the near-UV OH lines at 3060-3200 A are 
sensitive to the gas pressure and electron density in the 
model atmosphere. In contrast, the [Oi] line is roughly 
half as sensitive as the OH and 1 lines, and is insensitive 
to changes in the ratio of a-element to iron. Also, the OH 
lines are very crowded in the UV region, which makes their 
measurement difficult. 

Nissen et al. (|200C|) obtained the O abundances of 13 
dwarfs and subgiants using data from the ESO VLT. They 
suggested that the [O/Fe] is nearl y con stant at w +0.4 
for [Fe/H] below « -1.0. Maciel ( |2000[ ) suggested that 
[O/Fe] is not higher than +0.4 according to the radial 
O/H abundance gradients in Hn regions, hot stars and 
PNe, and [Fe/H] gradients from open cluster stars. 

From the foregoing discussion and considering the sim- 
ilar behavior between O and a elements ([a/Fe] is about 
-0.4 in me tal-po or stars, see Zhao & Magain 1990 and 



Fuhrmann 1998 ), we use the lower O abundances de- 
rived from the [Oi] line as observational constraints on 
our GCE calculations. And we will specially discuss the 
higher abundances in metal-poor stars obtained from the 
Oi triplet and OH lines in Sect. 8. Because the observed 
O abundances of metal-poor dwarfs are not enough to re- 
flect the abundance evolution, we have supplemented our 
data with some metal-poor giants. 

2.3. [N/Fe] vs. [Fe/H] 

The Ni line can be used in the abundance determina- 
tion of nitrogen, but this can only be done for reasonably 
metal-rich stars. Clegg et al. ( 1981 ) obtained the N abun- 
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dances of 20 disk stars (-0.9 < [Fe/H] < +0.4) using 
the Ni line, and suggested that [N/Fe] « 0.0. For halo 
stars, some results have been derived from NH lines in 
the UV region. Tomkin & Lambert fll984| ), Laird fl!985[ ) 
and Carbon et al. ( 1987] ) obtained similar results, [N/Fe] 
w 0.0, irrespective of metallicity in the range —2 < [Fc/H] 
< +0.3. Indeed, there are still some difficulties in deriving 
reliable N abundance from the NH bands at AA 3360 and 
3370 A. The difficulties are (a) uncertainties in the contin- 
uous opacities at these wavelengths, (b) severe atomic-line 
contamination of the NH features, (c) uncertainty in the 
dissociation energy of this molecule, (d) the crowdedness 
of the lines in the UV region. Also, it is very difficult to 
derive the N abundance from the CN bands at AA 4200 or 
3800 A. Possible reasons are that the dissociation energy 
of CN is not known to better than a factor of 2, and that 
the determination of N abundance from CN depends on 
a prior determination of reliable C abundance from other 
features, and that the CN bands are extremely weak in 
very metal-poor stars. A detailed discussions can be found 
in Wheeler et al. (|1989|). 



2.4- Table of observations 

Table 1 summarises the sources of the observations of 
carbon, nitrogen and oxygen. Most of the objects are 
dwarfs and subgiants; a few metal-poor giants provide ad- 
ditional O abundance data derived from the [Oi] A 6300 
A forbidden line. 



2.5. The age-metallicity relation and the metallicity 
distribution of G dwarfs 

The age-metallicity relation for the solar neighborhood is 
one of the important observational constraints; it shows 
the [Fe/H]-represented metallicity as a function of the 
star's age. The observational data and scatter are taken 
from Carlberg et al. ( 1985[ ), Meusinger et al. (1991) and 



Edvardsson et al. (1993) (their tables 14 and 15). 

Another important observational constraint is the 
metallicity distribution of long-lived G-type stars. Because 
these stars have long main-sequence lifetimes of about 15 
Gyr, comparable to the estimated age of the Galaxy, they 
represent a sample of stars which have not been removed 
by stellar evolution. The observational metallicity distri- 
bution of G dwarfs is taken from Chang et al. ( 1999[ ). 
The solar abundances of CNO elements are taken from 



Grevesse & Sauval (1998) 



3. Model of galactic chemical evolution 

A general, standard infall GCE model, namely, a multi- 
zone model for GCE without mass exchange among the 
zones is adopted in our present study. The rate of change 
of the surface mass density of isotope i in every zone in 



the gas at a given Galactocentric radius r and time t is 



dt 



0.5 



+ / X l (r,t-T M )^(r,t-T M )^(M)dM 

Mi 

M bu 

+ (3 J *(M 6 ) 

Mm 

M h „. 



f(p)Xi(r, t - T M2 )*(r, t - TM 2 )dfi 



dM h 



+ (1-/3) J X i (r,t-T M )^(r,t-T M )HM)dM 



+ / Xi(r,t-T M )y(r,t-T M MM)dM 



dai(r, t) 



dt 



(1) 



inf 



where the first term describes the disappearance of iso- 
tope i due to new star formation, the second term de- 
scribes the enrichment due to stellar wind mass loss of 
low mass single stars, the third term represents the rate 
of enrichment due to binary systems that undergo SN la 
explosions. The free parameter [3 is the amplitude factor 
of Type la supernova (hereafter SNIa) explosion and /(//) 
is the binary distribution function. The fourth term repre- 
sents the enrichment due to single stars or binary systems 
that do not undergo Type la events in the mass interval 
Mu to Mbu. The fifth term represents the enrichment rate 
due to massive stars that become SNII explosions with 
masses between Mbu to M u . The last term represents the 
infall of the primordial material. Xi(r,t) is the elemental 
yield from nucleosynthesis (see Sect. 4). More details have 
been given by Matteucci & Frangois ( 1989| ) and Timmes 
et al. ( |1995| ). 

• ^(r,t) (SFR) Analytical prescriptions for the SFR 
can be obtained from the observations. Schmidt (1959) 
suggested that the SFR is related to the surface mass den- 
sity of interstellar gas o~ g . Subsequent studies showed that 
the SFR depends not only on the surface mass density of 
gas but also on t he tot al surface mass d ensity of the ISM 
(Talbot & Arnett |l975[ Dopita & Ry der |l994[ Timmes et 
al. |1995| ; Chiappini et al. |l997| ; Dwek |l998|) . 

W e adopt a similar formula to that of Timmes et al. 
( |1995|) : 



*(r,t) = WtotM) 



.(r,t) 



o-tot (r, t) 



M pc- 2 Gyr- 



(2) 



where n is in the range 1-2 (Timmes et al. 1995), which 



we fix at 1.5 throughout our calculation, v is the efficiency 
factor in unit Gyr -1 , which can be adjusted around the 
value 1.0. 

• $(M) (IMF) For the solar neighborhood, the IMF 
is given in a s eries o f papers. The first formula was applied 
by Salpeter ( 1955 ), who suggested the single power-law 
form: $(M) ~ Af~( 1+X ) [x = 1.35), normalized to 1 over 
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Table 1. Observational data from different references and the line sources 



References C 



line, wavelength (A) 



N 

line, 



vavelength (A) 



O 

line, 



vavelength (A) 



[Fe/H] 



Numbers of stars 



1 


Ci, [Ci],CH Ni 8683 


Oi 6158, [Oi] 6300 


[-0.90 


+0.40] 


20 


2 


[C i] 8727 




[-1.00 


+0.25] 


85 


3 


[C fj 8727 




[-1.06 


+0.26] 


80 


4 


Ci 7100 




[-0.80 


+0.20] 


105 


5 


Ci 9100 


O i 7770 triplet 


[-3.00 


-0.80] 


34 


6 


CH 4300 NH 3360 




[-2.45 


+0.50] 


116 


7 


CH 4300 




[-2.60 


-0.70] 


32 


8 


CH 4300 NH 3360 




[-3.20 


-1.50] 


83 


9 


CH C 2 CN 


Oi, [Oi] 


[-2.61 


+0.12] 


19 (some giants) 


10 


NH 3360 




[-2.30 


-0.30] 


14 


11 




[O i] 6300 


[-1.60 


+0.30] 


7 


12 




[Oi] 6300 


-1.74, 


-1.43 


2 


13 




[O fj 6300 


[-0.80 


+0.30] 


23 


14 




[Oi] 6300 


-2.84, 


-2.31 


2 


15 




[O i] 6300 


[-1.80 


-0.70] 


13 


16 




O i 7770 triplet 


[-1.03 


+0.26] 


86 


17 




O i 7770 triplet 


[-1.00 


+0.10] 


90 


18 




O I 7770 triplet 


-1.84, 


-2.08 


2 


19 




O I 7770 triplet 


[-2.50 


-0.50] 


14 (some giants) 


20 




O I 7770 triplet 


[-2.72 


-0.47] 


14 


21 




O i 7770 triplet 


-3.00 




1 


22 




O i 7770 triplet 


[-3.50 


-0.20] 


30 


23 




OH 


[-3.00 


-0.30] 


24 


24 




OH 


[-3.02 


-0.50] 


24 


25 




OH 3138-3155 


[-3.20 


-1.80] 


9 


26 




OH , CO 


-1.22 




1 


27 




[O i] 6300 


[-1.20 


-0.10] 


24 (some giants) 


28 




[O f] 6300 


[-3.00 


-1.10] 


20 (giants) 


29 




[O i] 6300 


[-2.35 


+0.37] 


18 (giants) 


30 




[O fj 6300 


[-2.88 


-1.80] 


10 (giants) 



1. Cle gg et al. |1981| 2. Andersson 
Laird h.985|; 7. Tomkin et al 



Edvardsson 1994 



1991 ; 12. Spiesman fc Wallerstein 



Edva rdsson et al 
1993; ; 21. King |1994| ; 22. Abia 
Balachandran 
et al 



1991 



1987 



19861; 8. Carbon et al 

13. Nissen & Edvardsson 

18. Beveridge & Sneden 
23 



1991 



1993j; 17. Chen et al 

Rebolo 
Carney 1996]; 27. Barbuy 



200C 
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a given mass range. Then Scalo ( 1986 ) suggested that the 
power law should become steeper for M > IMq, x = 1.7. 
Following them, more authors suggested that IMF should 
flatten at the low mass end, for example Kroupa et al. 
( |l993|) . Dwek ( 1998 ) analyzed the se res ults, and gave a 
similar formula as Kroupa et al. ( 1993 ) (for three mass 
ranges). An important issue for the IMF is it being "top- 
heavy" , meaning that the formation frequency of massive 
stars in the early Galactic stage was higher than at the 
present. Recently, the question of a changing IMF as our 
Galaxy evolves has been examined. Chiappini et al. (2000) 
found that the metallicity distribution of G dwarfs is well 
explained by a constant IMF, rather than by a chang- 
ing one. Martinelli & Matteucci ( |2000 ) found that, with 
a changing IMF, even though the metallicity distribution 



of G dwarfs could be reproduced, other important obser- 
vations of [O/Fc] would then not fit so well. 



In our calculation, we adopt the formula of Dwek 
( |1998| ): 



$(M) 



C 3 M- 2(i 



for 0.1 < M/Mq < 0.3 
for 0.3 < M/Mq < 1 
for 1 < M/Mq < 100, 



(3) 



where C\, Ci and C3, the normalized coefficients, are de- 
Ma 

termined from the normalization / M$(M)dM = 1, 

Mi 

with Mi = O.IMq and M u = 100M Q . In fact, the 
bulk of chemical enrichment is caused by the stars with 
M > IMq. Thus, it is meaningful to fix the fraction £ 
of the total stellar mass distributed in stars above 1M 
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(Mi), which is equivalent to fixing Mi, namely, 

M u Mi M u 

M<P(M)dM = J M$(M)dM + J M$(M)dM 



M, 



Mi 
Mi 



M 1 



J M$(M)dM + ( = 1. 



(4) 



a = 2.50 is too high to fit the observations (M2K), so we 
only consider the results for a — 1.68, 2.00. 

4-2. Binary stars through SNIa explosions 

Type la SN explosions are assumed to occur in close bi- 



M, 



Pagel ( 19971) explained this method in detail. In our cal- 
culations for different sets of yields, we adjust ( so that 
the results give a good fit to the observations, especially 
the metallicity distribution of G dwarfs . 

•Infall The infall rate is generally taken to follow an 



exponent law (Matteucci & Frangois 1989; Timmes et al. 
1995|) or Gaussian form (Prantzos & Silk 1998) ; Chang et 



nary systems (Whelan & Iben 1973). In this model, the 
explosion is caused by a carbon-deflagration of the mate- 
rial accreting on the degenerate white dwarf (Nomoto et 
al. 1997b), and the ejecta is dominated by the 56 Fe iso- 
tope. In this paper, the nucleosynthesis yields of SN la ex- 
plosions are taken from the classical W7-model of Nomoto 
et al. dl997b|), in which ~ 0.613M Q of Fe is produced. 



al. |1999|) . We adopt the specific form of Timme s et al. 4-3- Massive single stars 



(|1995|) , but with o-q = 55M pc~ 2 (Sackett |l997| ) 



•Main-sequence lifetimes tm The main-se quenc e 
lifetimes of stars are taken from Schaller et al. ( 1992|) , 
which have been used generally in GCE calculations, such 
as Timmes et al. fll995| ) and Prantzos & Silk ( |1998[ ). 



4. Stellar yields and nucleosynthesis 

The stars are divided into binary and single stars. We 
only consider SNIa explosion for the binary stars, and 
divide the single stars into two groups: intermediate-, low- 
mass stars (0.9 — 8M , or 0.9 — 6M Q in case of strong 
overshooting) and massive stars (M > 8Af Q , or M > 6Mq 
in case of strong overshooting). 



4--1- Intermediate-, low-mass single stars 

The first detailed nucleosynthesis calculation for ILMS 
was made by RV (M = 1 — 8M Q , Z = 0.004, 0.02). Re- 
cently, VG extended the calculation to M = 0.8 — 8.0M Q 
and Z = 0.00 1, 0.004, 0.008, 0.02, 0.04. Forestini & Char- 
bonnel ( 1997| ) calculated the evolution of stars with M = 
3, 4, 5, 6, 7M Q and Z = 0.005, 0.02. MBC calculated 
the evolution of stars with Al = 0.819 — 5.0M Q and 
Z = 0.008, 0.02 for the case of strong overshooting, which 
means that the stars evolve rapidly so that the upper 
mass limit of ILMS is changed from the standard value 
of 8Mq to 6Mq. M2K gave a new set of stellar yields (for 
M = 0.817 - 5.0M@ and Z =0.004,0.008,0.02) based on 
the results of an updated evolution calculation. 

In our calculations, we choose the results of RV (their 
a = 1.5,7? = 0.333 case), VG, MBC and M2K as the 
yields of ILMS for their wide mass ranges and metallicity- 
dependent results. RV only gave results of 3.25M Q < M < 
8M Q for the case of a = 1.5, 77 = 0.333. In fact, the hot- 
bottom burning (HBB) process is not important for low 
mass stars, but is very significant for intermediate mass 
stars above 4M . So we choose the results with a = 0.0 
for the lower mass stars. M2K calculated three sets of 
yields with a = 1.68, 2.00, 2.50 for stars with M > 3.5M Q . 



A star in the mass range of 8— IOMq (or 6— 8M in the 
case of strong overshooting) generally develops a degener- 
ate O-Ne-Mg core after C-burning and eventually explodes 
as an electron capture supernova, leaving a neutron star of 
1.3Mq as a remnant and expelling a very small quantity 
of heavy elements. 

A star with a higher initial mass will end its life with 
an Fe-core collapse Type II SN explosion. Some authors, 
including WW, N97, M92 and PCB, calculated the stellar 
nucleosynthesis yields for such stars either with or with- 
out mass loss through wind. The results differ because 
of different choices of the parameters. We have carefully 
compared their parameters and results in Liang & Zhao 
( |200lD . Some of the parameters adopted by them are listed 
in Table 2, including the stellar mass, metallicity, and 
12 C(a,7) 16 reaction rate. 



Limongi et al. (2000) calculated the presupernova nu- 
cleosynthesis of massive stars with M — 13, 15, 20, 25M Q 
and Z = 0.02, and the explosive nucleosynthesis of stars 
with Z = 10~ 3 ,0. They did not give the yields of stars 
with M > 25Mq, which are very important for tracing 
the chemical evolution. And they did not give the explo- 
sive nucleosynthesis of stars with Z = 0.02. So we did 
not use their results in our calculation of the abundance 
evolution of CNO elements in the present paper. 

4.3.1. Yields calculated by WW 

WW calculated, for a large mass range, the metallicity- 
dependent yields of SN II explosion, and the results have 



been widely use d in G CE models (e.g. Timmes et al. |1995 
Chiappini et al. 1997 ) 

Timmes et al. (1995) suggested that the results of 



abundance evolution would be better if the Fe yields of 
WW is reduced by a factor of 2. Samland ( |1998[ ) sug- 
gested that the actual Fe yield of SNII explosion should 
be 0.046M Q , approximately a factor of 2 less than the WW 



value. Carigi fl994| ) used 0.075M Q as the Fe y ields o f SN II 
explosion in their GCE model. Chieffi et al. ( 1998 ) calcu- 
lated the nucleosynthesis of SNII explosion for a 25M Q 



8 Y. C. Liang et al.: Sources of carbon 

Table 2. Selected parameters of stellar evolution and nucleosynthesis of massive stars 



upper mass 
mctallicity 

stellar wind mass loss 
12 C(a,7) le O rate 



WW 

4OA/ 

Z = 1, 0.1, 0.01, 10~ 4 , 
No 

1.7xCF~0.74xCFHZ 



oz e 



N97 
70M o 

Zq 
No 

CFHZ 



M02 

12OM 

Z = 0.001, 0.02 

Yes 

CFHZ 



PCB 

12OA/ 

Z = 0.0004 - 

Yes 

CF 



0.05 



CFHZ=Caughlan et al. ( 138i 



CF=Caughlan & Fowler \ L98f ) 



star, and got an Fe yield of about 0.075M Q . SN 1987A, for 
an initial main sequence mass of 20M Q , ejected O.O75M 
of Fe, while for the same initial mass, WW got 0.151Af Q , 
again twice as large. Because of these findings we reduced 
the results of WW by a factor of 2 in our calculations. 

4.3.2. Yields given by N97 



N97 extended the calculations of Thielemann et al. (1996) 
to include 13 — 70M Q , solar metallicity stars. We should 
note that N97 used the evolution of the He core instead 
that of the entire star. They used a relation M(M a ) to 
transfer the result for a helium core of mass M a to that 
of a star of initial mass M (Sugimoto & Nomoto |1980p . 
This can be understood since the total ejecta of star must 
be the sum of the contributions of both the He-core and 
the envelope of mass M — M a . (Thomas et al. 1998). 



4.3.3. Yields given by M92 and PCB 

Both M92 and PCB considered radiative-driven stellar 
wind in their calculations, which was important for mas- 
sive stars, especially the W-R stars. The synthesized ele- 
ments of C, N and O, in particular, C, can be ejected into 
the ISM in stellar wind. 

How to choose the yields of W-R stars? Massey et 
al. (1995) studied the massive stars both in the general 



field and associations of the Magellanic Clouds. Their data 
showed that stars with initial masses M > 30M Q evolved 
through the W-R phase in the Large Magellanic Cloud 
(LMC); while the statistics of the Small Magellanic Cloud 
(SMC) are consistent with a somewhat higher mass limit 
of possibly 5OAf . Conti ( |1995| ) suggested that the W-R 
stars are highly evolved, luminous, hot, and (mostly) He- 
burning descendants of the most massive stars in the solar 
neighborhood, with masses M > 35M . In the present pa- 
per, we follow the suggestions of Maeder & Conti (1994) 
and the references therein and assume that the stars with 
M > 40Mq will evolve through the W-R stage. So we di- 
vide the massive stars into two mass ranges, M < 40M Q 
and M > 40A/q, in our discussion on the source of car- 
bon. For the upper mass limit of massive stars, we choose 
IOOMq, following the suggestion of Leitherer (1995), and 



the usual practice with the upper mass limit of IMF used 



in many GCE models (e.g. PCB, Prantzos & Silk |1998| ) 

M92 did not give Fe yields. When we use the M92 set 
of yields, we choose the following values: 0.15Mq for stars 



with initial main sequence mass M < 14M©, 0.075M Q for 
14M < M < 40M©, and 0.15M© for M > 40M Q . These 
choices are based on the following considerations. In the 
case of SN 1987A (2OAf during the main-sequence stage), 
the light curve, powered by the decay of 56 Ni and 56 Co, 
gives a determination of the produced 56 Fe of 0.075A/ Q , 
and for SN 1993J (main-sequence mass 14Af Q ), one of 
0.15A/ Q . (Thielemann et al. 1996). Carigi (1994) adopted 
O.O75M of Fe for a SNII explosion and 0.15Af Q for a 
SNIb explosion. Since SNIb takes place in W-R stars, 
and since we assume that the stars with M > 40Af© will 
undergo the W-R stage, so we use the value 0.15Af Q in 
such cases. 

M92 only gave the yield of N ejected in wind, and not 
the final value. Thus the N yields used in our calculations 
are lower limits. 

5. Results and analyses 

We set up a standard infall model of GCE to map out the 
abundance evolution of carbon, nitrogen and oxygen in 
the solar neighborhood. We choose the general formulae 
of SFR and IMF, and slightly adjust them by changing 
v and £ respectively. The 8 sets of nucleosynthesis yield, 
described in Sect. 4 are then used to set up 8 specific 
models. For each model, we try to fit the observational 
age-metallicity relation, the metallicity distribution of G 
dwarfs and the [O/Fe] vs. [Fe/H] relation. The calculated 
results will then be discussed: the [C/Fe] vs. [Fe/H] rela- 
tion, in this section, the [N/Fe] vs. [Fe/H] and [O/Fe] vs. 
[Fe/H] relations in Sect. 6 and Sect. 7, respectively. 

5.1. VG+WW 

In this model labelled (VG+WW), we use the yields of 
VG for ILMS, and WW for massive stars (see Sect. 4). 
The calculated age-metallicity relation, metallicity distri- 
bution of G dwarfs and [O/Fe] vs. [Fe/H] relation can all 
fit the observations well (Figs, la, c, d). The predicted 
abundance evolution of carbon is given in Fig. lb; this 
shows that [C/Fe] vs. [Fe/H] is approximately constant in 
time. The predicted trend agrees with the observations, in- 
cluding even the increasing positive [C/Fe] below [Fe/H] 
< —1.5. The reason of positive [C/Fe] may be that for the 
massive stars, the C yields are high and the Fe yields are 
relative low. In medium metal-poor region, [Fe/H] > —1.0, 
[C/Fe] slightly increases, showing the important C contri- 
bution from ILMS. Up to the late evolutionary stage of 
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the Galaxy, the ratio of SNIa rate to SNII rate increases, 
therefore the [C/Fe] value decreases with increasing [Fe/H] 
in the range of [Fe/H] > —0.5. The results show that the 
combined contribution of carbon by VG and WW can ex- 
plain the observed [C/Fe] vs. [Fe/H] relation. 

5.2. VG+N97 

Using the yields of N97 for massive stars and those of VG 
for ILMS (model VG+N97), our GCE results can fit the 
observational age-mctallicity relation, metallicity distribu- 
tion of G dwarfs and [O/Fe] values (Figs. 2a, c, d). But the 
predicted C abundance is too low to fit the observations 
(Fig. 2b). Compared to WW, the lower C yields calculated 
by N97 may be mainly caused by the choice of a higher 
rate of 12 C(a,7) 16 reaction. Also, the authors of N97 
themselves have compared the different nucleosynthesis 
results based on two differe nt 12 C(a,7) 16 reaction rates 
from Caughlan et al. (1985) (hereafter CFHZ) and Caugh- 
lan & Fowler fll988j ) (hereafter CF) respectively. CFHZ 
gave a higher reaction rate than CF: CFHZft:2.3xCF; 
and WW adopted 1.7xCF«0.74xCFHZ in their calcu- 
lation, which is lower than in N97. So when [O/Fe] can 
be matched, [C/Fe] is lower than the observations. These 
results underline the importance of the 12 C(a,7) 16 re- 
action rate to elemental nucleosynthesis. 

Certainly, besides the 12 C(a,7) 16 reaction rate, con- 
vective mechanism and explosion energy can affect the 
nucleosynthesis yields (N97). 



5.3. VG+M92 

The results predicted by model VG+M92 can fit the ob- 
served age-metallicity relation, metallicity distribution of 
G dwarfs and [O/Fe] vs. [Fe/H] (Figs. 3a, c, d). There is a 
strong increase in the calculated [C/Fe] for [Fe/H] > — 1.2 
(Fig. 3b). What is the main contributor? ILMS? Stellar 
winds of W-R stars? Or a combined contribution of the 
two? Let us try to identify the main source. 

Fig. 3e gives the result for only the yields of massive 
stars of M92. It shows the predicted [C/Fe] increases with 
metallicity in metal-rich region. The calculated [C/Fe] 
using the yields of ILMS (VG) and only massive stars 
with M < 40M Q (M92) is presented in Fig. 3f: here, the 
calculated [C/Fe] is slightly lower than the observations 
in metal-rich region. The difference between Fig. 3e and 
Fig. 3f shows that the contribution to carbon from stellar 
wind of high metallicity W-R stars calculated by M92 is 
greater than the contribution from ILMS given by VG. 
The M > 40Mq massive stars with higher metallicity, 
the W-R stars, can eject significant amounts of carbon, 
which causes the [C/Fe] value to increase with increasing 
metallicity in the range of [Fe/H] > —1.0 in Fig. 3e. 



5.4. RV+WW 

Using results of RV as the yields of ILMS and WW for 
massive stars (model RV+WW), we calculated the cor- 
responding age-metallicity relation, metallicity distribu- 
tion of G dwarfs, [O/Fe] vs. [Fe/H] and [C/Fe] vs. [Fe/H] 
(Figs. 4a, c, d, b). When these three results match the 
observations, [C/Fe] first increases slightly and then de- 
creases, and is slightly higher than the prediction of model 
VG+WW in metal-rich region. These results show that 
the RV calculation gives more carbon than does VG due 
to the use of different parameter values, (see M2K, Liang 
& Zhao [2001| for details). 



5.5. RV+N97 

Figs. 5a, c, d show that the calculated results using model 
RV+N97 can fit the observations, i.e., the age-metallicity 
relation, the metallicity distribution of G dwarfs and 
the [O/Fe] vs. [Fe/H] relation. The calculated [C/Fe] vs. 
[Fe/H] then shows that more carbon is given by RV than 
VG, and this increases the predicted [C/Fe] in metal- rich 
region (Fig. 5b). Because N97 gives low carbon yields for 
the reason of adopting a low 12 C(a,7) 16 reaction rate, 
[C/Fe] falls below the observations in the whole metallicity 
range (Fig. 5b). 

5.6. RV+M92 

Figs. 6a— d display the results based on the yields of RV 
for ILMS and of M92 for massive stars (model RV+M92). 
Figs. 6a, c, d show the calculations can fit the correspond- 
ing observations, while Fig. 6b shows that [C/Fe] increases 
with metallicity in metal-rich region (Fig. 6b). There are 
no apparent difference here between this model and the 
model VG+M92 (Fig. 3b); this can be understood be- 
cause the C contribution from W-R stars given by M92 is 
higher than that from ILMS given by RV or VG. 

Fig. 6e shows the [C/Fe] vs. [Fe/H] relation based only 
on the yields of massive stars from M92; with no contribu- 
tion from ILMS, this is still higher than the observations 
in metal-rich region. Fig. 6f exhibits the predicted [C/Fe] 
using the yields of RV for ILMS and those of M92 for 
massive stars with M < 40M Q ; it fits the observations, 
and shows a slightly higher [C /Fe] than is given by model 
VG+M92 in metal-rich region (Fig. 3b). This result shows 
that RV gives higher C yields than does VG with a differ- 
ent choice of the parameters. 

5.7. MBC+PCB 

It is proper to combine the yields of MBC and PCB 
together since they used the same parameters (Padova 
group). They adopted a larger convective overshooting, so 
that stars with M > 6M© rather than the standard >8M Q 
would end their lives through SNII explosions. Both the 
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Fig. 1. Results using the nucleosynthesis yields of VG for ILMS and WW for massive stars (the solid lines) (model 
VG+WW). (a) The agc-metallicity relation in the solar neighborhood; (b) The [C/Fe] vs. [Fe/H] relation; (c) The 
metallicity distribution of G dwarfs in the solar neighborhood; (d) The [O/Fe] vs. [Fe/H] relation. The corresponding 
observations are shown in the figures. 



age-metallicity relation and metallicity distribution of G 
dwarfs can fit the observations well (Figs. 7a, c). The O 
abundance is slightly lower than the observations (Fig. 
7d), which may be caused by the lower 12 C(a, 7) 16 reac- 
tion rate (CF) used by the Padova group. The predicted 
[C/Fe] vs. [Fe/H] relation is given in Fig. 7b, which shows 
a strong increase of [C/Fe] with increasing metallicity from 
[Fe/H] » -0.7 on. 

Fig. 7e displays the [C/Fe] based only on the contribu- 
tion of carbon from massive stars; it shows that W-R stars 
with a high metallicity (Z = 0.008, 0.02) can eject signif- 
icant amounts of carbon into the ISM, much more than 
the lower metallicity stars. Thus, [C/Fe] increases steeply 
starting from [Fe/H] w —0.7. The subsequent decrease 
from [Fc/H] « —0.3 on is the result of Fe contribution 
from SNIa explosions. 

Fig. 7f shows the case of only ILMS and the M < 
40M Q part of massive stars; there is basic fit to the obser- 
vations. Comparing Figs. 7b, 7e and 7f, we find the main 
effect of the C yields from metal-rich W-R stars is the ob- 



vious increase of [C/Fe] in metal- rich region shown in Fig. 
7b; — an increase to above the observations. 

5.8. M2K+PCB 

We use the nucleosynthesis yields of M2K in our GCE cal- 
culations because this new set of yields include the results 
of ILMS with the low metallicity value Z = 0.004, which 
may contain significant information on the ILMS during 
the early Galactic stage. 

We combine M2K and PCB in our model M2K+PCB. 
The calculated age-metallicity relation, metallicity distri- 
bution of G dwarfs and [O/Fe] vs. [Fe/H] relation again 
match the observations, and are similar to the results of 
model MBC+PCB, and so are not shown here. The pre- 
dicted evolution of [C/Fe] with metallicity is given in Fig. 
8a, which shows a strong increase with increasing metal- 
licity from [Fe/H] w —1.5 on. In order to understand this 
increase, we re-calculate the [C/Fe] vs. [Fe/H] relation us- 
ing only the contribution of massive stars (Fig. 8c), and 
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Fig. 2. Same as Fig. 1, but with the nucleosynthesis yields of VG for ILMS and N97 for massive stars (model VG+N97). 



only the yields of ILMS and M < 40M Q massive stars 
(Fig. 8d) respectively. Comparing Figs. 8c and 7e, and 
Figs. 8d and 7f, we guess that maybe it is the ILMS with 
Z = 0.004 that produce a significant amount of carbon, 
which lead to the strong increase of [C/Fe] displayed in 
Fig. 8a. Fig. 8e confirms our guess, in which the solid line 
represents the calculated [C/Fe] vs. [Fe/H] relation using 
only the yields of Z = 0.008,0.02 (the same metallici- 
ties as in MBC) of M2K for ILMS and PBC for massive 
stars, and the dashed line represents the results of model 
MBC+PCB (the same line as in Fig. 7b). 

The solid and dashed lines in Fig. 8f represent the cor- 
responding results for nitrogen. Figs. 8e, f show that the 
nitrogen yields of M2K with Z = 0.008, 0.02 are very simi- 
lar to the corresponding results of MBC, though M2K gave 
more carbon and nitrogen according to the updated pa- 
rameters, resulting in the slightly higher [C/Fe] and [N/Fe] 
values in metal-rich region. Fig. 8b illustrates the [N/Fe] 
vs. [Fe/H] relation using the yields of M2K for ILMS; the 
relation is higher than that shown in Fig. 8f. 

These results show that ILMS with low metallicity Z = 
0.004 contribute very important amounts of C and N to 
the ISM (also see Fig. 1 of M2K). The lower metallicity 



ILMS have higher C and N abundances because for them 
the HBB process is more efficient, and mass loss is less so. 
The lower mass loss rates correspond to longer TP-AGB 
lifetimes, hence a greater number of the third dredge-up 
events and a longer duration of the HBB process. The 
up-to-date nucleosynthesis calculations (VG) only go as 
low as Z = 0.001. Calculations of ILMS of even lower 
mctallicities are needed; these may revise the trends of 
elemental abundance evolution, and provide different hints 
to the early GCE history. 

Though the nucleosynthesis results with a = 1.68 and 
2.00 given by M2K are different (more nitrogen is pro- 
duced by intermediate mass stars with a — 2.00 than with 
a = 1.68), we find no obvious differences between these 
two sets of results. So we only give the GCE results based 
on the yields with a — 2.00 in this paper. 

In summary, we give and analyze the calculated abun- 
dance evolution of carbon using 8 recently published sets 
of nucleosynthesis yields. The parameters of the 8 models 
are listed in Table 3. Here are some of the obvious im- 
plications of the parameters. The higher the £ value, the 
more contribution comes from massive stars, and then the 
predicted O abundance will be higher, especially in low 
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Fig. 3. (a), (b), (c) and (d) are same as in Fig. 1, but with the nucleosynthesis yields of VG for ILMS and M92 for 
massive stars (model VG+M92). (e) The calculated [C/Fe] vs. [Fe/H] only considering the C yields of massive stars 
of M92; (f) The calculated [C/Fe] vs. [Fe/H] considering C yields of VG for ILMS and C yields of only M < 4OM 
stars from M92 for massive stars. 

mctallicity region. A higher (3 value means a higher SNIa gion. At the same time, a higher v will lead to a higher 
rate, hence a greater amount of Fe element from SNIa ex- peak in the mctallicity distribution of G dwarfs, 
plosion. A higher v value represents a higher SFR, more 
stars form, leading to an increase in the O abundance and 
[Fc/H], and the effect is more obvious in metal-poor re- 
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O Edvardsson et al. 1993 
^ Meusinger et al. 1991 
A Carlberq et al. 1985 



time (Gyr) 



(b) 




- 1 
[Fe/H] 



Chan" et at. 1999 



0.05 - 




(d) - 




-0.5 
[Fe/H] 



-1 
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Fig. 4. Same as Fig. 1, but with the nucleosynthesis yields of RV for ILMS and WW for massive stars (model 
RV+WW). 

Table 3. Input parameters of our 8 models 



model 


c 


P 


V 


model 


c 


9 




VG+WW 


0.51 


0.037 


1.0 


RV+WW 


0.54 


0.033 


1.0 


VG+N97 


0.40 


0.045 


1.0 


RV+N97 


0.40 


0.040 


1.0 


VG+M92 


0.40 


0.040 


0.7 


RV+M92 


0.40 


0.035 


0.7 


MBC+PCB 


0.38 


0.040 


0.5 


M2K+PCB 


0.38 


0.040 


0.5 



5.9. The 12 C(a,7) 16 O reaction rate 



The 12 C (a, 7) 16 reaction rate is very important to stellar 
nucleosynthesis. It determines the C and O abundances, 
hence the abundances of the heavier elements. Generally, 
in nucleosynthesis calculations, the 12 C(a,7) 16 reaction 
rates given by CFHZ and CF are widely used. CFHZ gives 
higher values than does CF: CFHZ«2.3xCF. And the 
higher reaction rate leads to a lower C abundance. N97 
exhibited this result in their Fig. 7. They suggested that 
the actual reaction rate was between CFHZ and CF. WW 
adopted 1.7xCF («0.74xCFHZ) to be the actual reac- 
tion rate in their nucleosynthesis calculations. They gave 
higher C yields than did N97. 



For the higher massive stars, stellar wind mass loss 
strongly affects the C yields, so the effect caused by the 
12 C (a, 7) 16 reaction rate is obscured. The lower reaction 
rate of CF has led to the lower O yields of PCB. 

5.10. Understanding the W-R stars 



Leitherer (1995) suggested that the most massive stars 
known in the universe have masses around 100M© and 
lifetimes of a few Myrs. They are rare: in the solar neigh- 
borhood only about one such massive star is counted per 
10 5 to 10 6 solar-type stars. There are three key phases 
during the evolution of a massive star: O star, luminous 
blue variable star and Wolf-Rayet star. The mass loss of 
a massive star can be calculated from UV wind lines, H a 
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Fig. 5. Same as Fig. 1, but with the nucleosynthesis yields of RV for ILMS and N97 for massive stars (model RV+N97). 



or radio fluxes. The UV P Cygni profile, ubiquitous in 
O-type stars, provides a direct indication of stellar wind; 
H a has been recognized as the prime source of mass-loss 
in early-type stars; and winds in hot stars can be read- 
ily observed at IR-mm-radio wavelengths via the free-free 
"thermal" excess caused by the stellar wind. The mass 
loss rate scales as Z 5 for O stars and early B stars, as 
Z 8 for mid-B sup ergia nts, and as Z fo r A supergiants 
(Kudritzki & Puis fcOOOt Crowthcr |2000a|) . 

W-R stars are easily detectable due to their intrinsic 
brightness and the profound effects on the interstellar en- 
vironment. They eject a large amount of material into the 
ISM via the radiative-driven stellar wind. The nucleosyn- 
thesis calculations show that the ejected carbon by mas- 
sive stars are very important to the enrichment of the ISM. 
In general, the mass loss rate of massive stars adopted 
in the stellar evolution calculations are M ~ M 2 5 and 
M ~ Z 5 (M92; PCB). However, according to recent re- 
views, our knowledge of the mass loss of massive stars is 
not very clear. 

Meynet & Maeder (2000) calculated the evolution of 



and N— enhancements; the enhancements are the greater 
for the hig her th e mass and the rotation. In his review, 
Crowther ( 2000a ) found that the recent mass- loss rates 
for Galactic W-R stars indicated a downward revis ion of 
2—4 relative to previous calibrations (e.g. Langer 1989 : 
M ~ M 2 - 5 ) due to clumping (e.g. Schmutz 1997|) . Indeed, 
metallicity affects the yields of carbon. When mass loss 
gets very efficient, 12 C will increase because He-burning 
products can be revealed on the surface. If mass loss 
is extreme, anyway, 12 C and 16 yields may even de- 
crease because most of the mass is rapidly lost in the 
wind in the form of 4 He. But at present, the question 
of a Z— dependent mass loss in W-R stars remains open. 
Crowther ( 2000a ) suggested that there is a clear distinc- 
tion among spectral subtypes in galaxies with different 
metallicities, which is qualitatively explained by evolu- 
tionary and spectroscopic models (Crowther 2000b| ). WN 
stars in LMC show negligible spectroscopic difference from 
their Galactic counterparts. And the situation is less clear 
for the SMC, since most W-R stars are complicated by 
binarity (Crowther 2000a). 



stars with rotation, and suggested that rotation can in- 
crease mass loss, also can bring significant surface He— 
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Fig. 6. Same as Fig. 3, but with nucleosynthesis yields of RV for ILMS and M92 for massive stars (model RV+M92). 



5. 11. Contributions by stars of different masses 



Comparing our calculated results using different nucle- 
osynthesis yields (Sect. 5.1—5.8), considering the uncer- 
tainty of stellar wind mass loss of massive stars (see Sect. 
5.10), and assuming that the actual 12 C(a,7) 16 reaction 
rate to be between the values of CFHZ and CF, we choose 



to use the results of model VG+WW in an analysis of the 
contributions of carbon by stars of different masses. 

Here, the reason of choosing VG rather than RV for 
ILMS is because of suggestions that the yields of RV are 
not successful in predicting the observed C/O and N/O 
ratios in stars and planetary nebulae, or in reproducing 
the carbon stars of relatively low lumin osities in the Mag- 
ellanic Clouds (Gustafsson et al. 1999). 
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Fig. 7. Same as Fig. 3, but with nucleosynthesis yields of MBC for ILMS and PCB for massive stars (model 
MBC+PCB). 



Fig. 9 shows the fraction of carbon contributed by mas- 
sive stars based on the WW yields for massive stars and 
the VG yields for ILMS, as a function of the Galactic age. 
In the early Galactic stage, almost all the carbon is pro- 
duced by massive stars. As the age advances, ILMS con- 
tribute more and more, and an overtake the massive stars 
from about 1.65 Gyr on. At the present epoch, many low 



mass stars have also evolved and ejected matter into the 
ISM, and the relative contribution by the ILMS is further 
increased. 
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time (Gyr) 

Fig. 9. The relative contribution of massive stars to car- 
bon production at different Galactic ages using the results 
of WW for nucleosynthesis yields of massive stars and VG 
for ILMS. 



6. Nitrogen Abundance 

The element nitrogen includes a primary and a secondary 
component. If nitrogen production depends on the ini- 
tial CNO abundance, the [N/Fe] would be dependent on 
[Fe/H], and this part of nitrogen is secondary; however, if 
nitrogen production is chiefly primary, then [N/Fe] would 
be independent of the value of [Fe/H]. 

Laird fll985| ) and Carbon et al. fll987p found the ob- 
served [N/Fe] to be constant over the entire metallic- 
ity range, although the precise level is uncertain. So, it 
seems that the observations require a primary nitrogen 
source from massive stars. However, if the observed con- 
stant nitrogen in low metallicity region is correct, then the 
present yields by massive stars cannot s upply the needed 
prim ary n itrogen source (Timmes et al. 1995 Chiappini 
et al. 1997 ). Our results, based on the recently published 
nucleosynthesis yields of massive stars, also point to an ab- 
sence of a primary nitrogen source. Figs. 10a, b, c show 
the calculated results of [N/Fe] vs. [Fe/H] using our GCE 
models VG+WW, VG+N97 and VG+M92 respectively. 
Those models with the RV yields of ILMS give similar 
trends to these with the VG yields; so we do not show 
them here. The [N/Fe] values of model MBC+PCB and 
M2K+PCB will be given in Sect. 5.8. 

In these figures, the solid lines represent the results 
with the N yields of both ILMS and massive stars, and the 
dashed lines represent the results with only the contribu- 
tion from massive stars. These results show that ILMS are 
the dominant donor of nitrogen, while massive stars have 
provided only a small part to the ISM during the whole 
Galactic age up to the present. However, in the range of 
—3 < [Fc/H] < —2, [N/Fe] increases steeply (solid lines in 
Figs. 10a, c), due to intermediate mass stars as source of 
secondary nitrogen. As metallicity increases up to [Fe/H] 



f» —1.0, [N/Fe] increases gradually, — typical behaviour of 
secondary nitrogen. As [Fc/H] increases further, [N/Fe] 
approaches the solar value, which shows that the primary 
nitrogen from LIMS plays an important role in the pro- 
duction of nitrogen since that time. 

From another aspect, comparison between the results 
shown in Fig. 10a and Fig. 10c confirms further that the 
ILMS are the main source of nitrogen. These two figures 
are based the models VG+WW and VG+M92 which use 
different N yields from massive stars, yet the results are 
very similar, so we infer that it is the N yields of ILMS 
that provide the bulk of nitrogen to the ISM. In addi- 
tion, we should note that M92 gives only the N yields of 
stellar wind and does not include contribution from SN II 
explosions, so the dashed line is higher in Fig. 10a than 
in Fig. 10c. Again, the predicted N abundance evolution 
by the models MBC+PCB and M2K+PCB also support 
the statement that ILMS are the dominant N contributor 
(see Figs. 8b, f). The model VG+N97 predicts high [N/Fe] 
in early Galactic stage, though not as high as the solar 
abundance. A possible reason is that N97 only gives the 
yields of solar metallicity massive stars, which eject more 
nitrogen than do low metallicity stars. With increasing 
metallicity, Fe from SN la explosions plays an increasingly 
more important role in the enrichment of the ISM, and 
this leads to a [N/Fe] decrease with the Galactic evolu- 
tion (dashed line in Fig. 10b). 

In general, the relative weight of the secondary and pri- 
mary components on the theoretical yield of 14 N of ILMS 
depends on the interplay between the secondary enrich- 
ment caused by the first and second dredge-up episodes, 
and the primary contribution given by the CNO-cycle dur- 
ing the envelope burning. 13 C and 14 N are regarded as sec- 
ondary elements in the sense that they are formed from 
12 C and 16 O originally present in the star at the time of its 
formation. Typical in this sense is the secondary produc- 
tion of 13 C and 14 N during the first dredge-up, and the 
secondary production of 14 N during the second dredge- 
up. However, when the third dredge-up process operates 
in conjunction with the HBB process, primary 13 C and 
14 N are produced. Meanwhile the primary 12 C, dredged 
up after each He-shell flash is later converted in part into 
primary 13 C and 14 N by the burning at the base of the 
convective envelope during the interpulse phase. 

Is it sure that the observed [N/Fe] is near 0.0 in 
metal-poor region? Carbon et al. (1987) found two stars 
with high N abundances, HD 74000 with [N/F e] - +0.9, 

HD 25329 with [N/Fe] h0.45. Laird ( |1985| ) obtained 

high N abundances for four stars, HD 74000 with [N/Fe] ~ 
+1.15 and [Fe/H] = -2.02, HD 97916 with [N/Fe] - +0.75 
and [Fe/H] = -1.30, HD 160617 with [N/Fe] - +1.65 and 
[Fe/H] = -1.99, HD 166913 with [N/Fe] - +0.70 and 
[Fe/H] = —1.80. The possible reasons for the very high 
N abundance given by the authors are: (1) internal mix- 
ing; (2) mass transfer in a binary system; (3) primordial N 
enhancement. Beveridge & Sneden (1994) reanalyzed the 



18 



Y. C. Liang et al.: Sources of carbon 



chemical abundances of two N-rich dwarf stars, namely 
HD 74000 and HD 25329. Their results showed that all the 
very heavy elements synthesized through neutron-capture 
s-process were enhanced in these two stars. It is likely 
that the enrichment of N and the very heavy elements 
originated in material dredged up from the helium burn- 



[O i] forbidden line (+0.4) are lower than those of dwarfs 
from the Oi triplet at A 7770 A and OH lines (+1.0) (see 
Fig. 11). And using the A 6300 A forbidden line, similar 
abundance values to the giants have been obtained for 
some metal-poor dwarfs, though the line is very weak in 



ing shell of former AGB stars. Also, Carbon et al. (1987) 
pointed out that, omitting the N-rich objects, the observed 
[N/Fe] vs. [Fe/H] plot showed a negative correlation with 
[Fe/H]. In this sense, our results can basically fit the ob- 
servations. More explorations of metal-poor N-rich objects 
are needed. 

We now follow Sect. 5.11 and calculate the contri- 
bution by massive stars to nitrogen using the model 
VG+WW. Fig. lOd shows the fractional contribution as a 
function of the Galactic age. We note that, except at the 
extreme early stage before the intermediate mass stars 
began to evolve, the relative contribution from massive 
stars is very low in the early stage, and that it gradually 
increases with the evolution of the Galaxy, then slightly 
decreases in the late stage. This trend means that the yield 
of nitrogen varies strongly with the metallicity. However, 
the proportion of nitrogen contributed by massive stars is 
generally very low, even the peak is less than 20%. WW 
calculated a very low N yield for stars with zero metal- 
licity, on the order of 10~ 6 M Q . The N yields of ILMS 
with low metallicities are higher than WW, for example, 
the newly formed and ejected nitrogen from 8M@ stars 
with Z = 0.001 is about 1O _2 M , thus, the intermedi- 
ate mass stars have been making important contributions 
since they began to evolve and eject materials into the 
ISM. With increase of metallicity, the N yield of massive 
stars also increases, as is shown by rising part of the curve 
of Fig. lOd. The eventual decrease in the curve at the late 
stage is probably caused by the higher N yield from the 
ILMS with higher metallicities. 



7. Oxygen Abundance 

The source of oxygen is generally clear, namely, it is pro- 
duced mainly in short-lived, massive stars and ejected 
through SNII explosions. As our Galaxy evolves, more 
and more Fe is produced through long lifetime SN la explo- 
sions, and the ratio [O/Fe] decreases. This trend is shown 
in all of our calculated results given in Sect. 5. 



Conti et al. (1967) presented the first data indicat- 
ing the existence of an oxygen overabundance in metal- 
poor stars. Many authors continued to explore the evolu- 
tion of [O/Fe], and confirmed the overabundance of oxy- 
gen in metal-poor stars (see Sect. 2). Edvardsson et al. 
( |1993| ) and Chen et al. ( 2000D separately obtained the 
oxygen abundances of a large sample of disk stars, both 
showed that [O/Fe] decreases with increasing metallicity 
in Galactic disk stars. For halo stars, however, available 
observations exhibit an apparent contradiction, namely, 
the [O/Fe] values of metal-poor giants derived from the 



these stars (Nissen et al. 2000). 

The difference can be understood from two aspects. 
In the first aspect, we believe that oxygen has a similar 
behavior to other a elements (Mg, Si, S, Ca), that is, we 
believe that the lower [O /Fe] in metal-poor region is true, 
and the relatively higher ratio derived from the O I triplet 
or OH line is wrong. The reasons are as follows. The Oi 
triplet has a very high excitation feature (9.15 eV) and is 
quite sensitive to the temperature structure of the adopted 
model atmosphere, and NLTE effects cause the higher O 
abundance derived from it. And the OH lines in the UV 
are very sensitive to the temperature structure of the stel- 
lar atmosphere and/or to NLTE effects too. According to 
this argument, there should be no obvious change in the 
O abundance as the dwarfs evolve to the giants. 

On the other hand, if we believe that the higher [O /Fe] 
derived from the O I triplet and OH line reflects the true 
O abundance, then we have to reconsider the abundance 
evolution of oxygen, hence the GCE result. Maybe the top- 
heavy IMF is true, which favours the formation of massive 
stars in the early Galactic stage, hence a high O abun- 
dance. This assumption would imply a corresponding high 
[a/Fe] in metal-poor dwarfs, which is not found in present 
available data. Could the difference between [O/Fe] and 
[a/Fe] mean that their behaviors are not similar? Per- 
haps NLTE effects can revise the present [a/Fe] value. 
Recently, Zhao & Gehren ( [2000D reported that NLTE ef- 
fects can increase the [Mg/Fe] abundance, and the correc- 
tions are larger for the lower metallicities of the metal-poor 
stars. Thus, maybe NLTE corrections can also increase 
the abundances of other a elements at the early Galactic 
stage. 

Goswami & Prantzos (2000) pointed out three reasons 
for the higher [O/Fe] value in low metallicity region: (1) 
The Fe produced by SN la explosions entered the galactic 
scene as early as [Fe/H] ~ —3, instead of [Fe/H] ~ — 1 
in the "standard" scenario; (2) the O yields from mas- 
sive stars are, for some reasons, metallicity dependent; (3) 
there is the possibility of the yields of Fe and all elements 
heavier than oxygen being metallicity-dependent. 

8. Conclusions 

Using the standard GCE model, we have attempted an 
exploration of the questions of the sources of carbon and 
the abundance evolution of CNO elements. Because the 
yield of stellar nucleosynthesis is a very important pa- 
rameter of the GCE model, and different yields give dif- 
ferent evolutionary behaviors of the elements, we chose 
the recently published different nucleosynthesis yields: the 
yields given by RV, VG, MBC and M2K for intermediate- 
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[Fe/H] 

Fig. 11. The observed [O/Fe] vs. [Fe/H] data. Oxygen 
abundances derived from [O I] 6300 A forbidden line are 
indicated by triangles, from O I triplet by pluses, and from 
OH lines in UV region by diamonds. 



low-mass single stars; the results of the classical W7 



model of Nomoto et al. (1997b) for binary stars through 
SNIa explosions; the yields given by WW, N97, M92 and 
PCB for massive single stars. Then we set up 8 spe- 
cific models for the combinations, VG+WW, VG+N97, 
VG+M92, RV+WW, RV+N97, RV+M92, MBC+PCB 
and M2K+PCB. 

Our results show that, after fitting the predicted age- 
metallicity relation, metallicity distribution of G dwarfs 
and [O/Fe] vs. [Fe/H] relation to the observations, very 
different [C/Fe] vs. [Fe/H] relations are given by the dif- 
ferent models with different yield combinations. 

The 12 C(a,7) 16 reaction rate is very important to 
stellar nucleosynthesis calculations. Hence, it affects the 
results of the GCE calculations. A higher rate of this re- 
action will produce less carbon, and a lower rate, more. A 
suitable choice is between the values given by CFHZ and 
CF. 

Generally speaking, in the early stage of our Galaxy, 
massive stars are the main contributor of carbon. As our 
Galaxy evolves to the late stage, the long lifetime ILMS 
begins to play an important role in the enrichment of 
the ISM. At the same time, the high metallicity W-R 
stars eject significant amounts of carbon into the ISM by 
radiative-driven stellar wind; the metallicity influences the 
outer opacities and therefore the mass loss rate by stellar 
winds in massive stars. However, on the present published 
nucleosynthesis yields, our results still cannot distinguish 
the source of carbon in metal-rich Galactic stage. It may 
be that just the massive stars (M > 8A1©) are sufficient to 
account for the observed [C/Fe] values (see Figs. 3e, 6e); 
or it may be just the ILMS and M < 40M© massive stars 
without the W-R stars are enough to explain the observed 
carbon in metal-rich stars (see Figs, lb, 7f). 



However, we can at least distinguish the different con- 
tributions of ILMS and massive stars to carbon from our 
calculations. Having comparing the GCE results shown 
in Sect. 5.1—5.8, and considered the undetermined char- 
acters of stellar wind mass loss of massive stars, and as- 
suming the actual 12 C(a,7) 16 reaction rate to be be- 
tween CFHZ and CF, we calculated the respective con- 
tributions to carbon from massive stars and ILMS in the 
model WW+VG. Fig. 9 shows that, in the early stage, al- 
most all carbon is produced by massive stars; with increase 
of the Galactic age, the contribution from ILMS becomes 
more and more and eventually exceeds that of the mas- 
sive stars. The latter contribution is increased when the 
contribution from W-R stars is included. 

Henry et al. fl2000] ) gave an excellent analysis on the 
source of C and N. But they only gave detailed results for 
the yields of VG for ILMS and M92 for massive stars. And 
they increased the yields of M92 in their calculations (see 
their Table 3), which increased the contribution of carbon 
from massive stars. 



Carigi ( 2000 ) compared different set of yield, and sug- 



gested that massive stars with stellar wind are the main 
source of carbon. Compared to their work, we use the 
abundances of a large sample of dwarf stars as one of the 
most important observational constraints rather than just 
the stellar age. The reason is that few halo stars have de- 
termined ages, while the results reported by Edvardsson et 
al. ( |1993| ), Gustafsson et al. ( |l999D and Chen et al. ( |2000[ ) 
are only for disk stars, so the element abundance in the 
early stage of the Galaxy cannot be tested by observations. 
In addition, we divided the massive stars into two ranges, 
Al < 40M Q and M > 40M o , so that we can compare 
the carbon contributions given by W-R stars and ILMS 
more clearly. And we use the new set of nucleosynthesis 
yields of M2K, which contain important information on 
low metallicity ILMS. 

Our results confirm that most of nitrogen is produced 
by ILMS, especially by intermediate mass stars. And the 
present yields given by massive stars cannot supply the 
needed primary N source. Massive stars produce the main 
part of oxygen in the universe. But the observed O abun- 
dances derived from different atomic lines lead to obvi- 
ous differences in the [O/Fe] determination in metal-poor 
stars. Maybe we should check the procedure and many 
corresponding parameters to obtain abundance of oxygen. 
Also, perhaps this difference hints at some new informa- 
tion on Galactic chemical evolution, especially, it is nec- 
essary to understand the evolution of our Galaxy in the 
early stage carefully. 
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Fig. 8. (a) The calculated [C/Fe] vs. [Fe/H] using nucleosynthesis yields of M2K for ILMS and PCB for massive stars 
(model M2K+PCB); (b) Same as (a) but for nitrogen; (c) The calculated [C/Fe] vs. [Fe/H] considering only yields of 
massive stars from PCB without ILMS; (d) The calculated [C/Fe] vs. [Fe/H] considering yields of M2K for ILMS and 
only M < AQMq massive stars from PCB; (e) The solid line represents the [C/Fe] vs. [Fe/H] relation considering the 
yields of PCB for massive stars and M2K only with Z = 0.008 and Z = 0.02 metallicities for ILMS, and the dashed 
line represents the same result as in Fig. 7b; (f) Same as (e) but for nitrogen. 
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Fig. 10. (a) The calculated [N/Fe] vs. [Fe/H] relation using nucleosynthesis yields of VG for ILMS and WW for massive 
stars. The solid line represents the result with both the contributions of ILMS and massive stars, and the dashed line 
represents the result with only massive stars, (b) Same as (a) but with the yields of N97 for massive stars, (c) Same as 
(a) but with the yields of M92 for massive stars, (d) The relative contribution of massive stars to nitrogen production 
with the Galactic ages predicted by using the yields of VG for ILMS and WW for massive stars. 



